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ABSTRACT 

We present new calibrations of the near-infrared surface brightness fluctuation (SBF) distance method for 
the El low (7iio) and F160W (//leo) bandpasses of the Wide Field Camera 3 Infrared Channel (WFC3/IR) on 
the Hubble Space Telescope. The calibrations are based on data for 16 early-type galaxies in the Virgo and 
Fornax clusters observed with WFC3/IR and are provided as functions of both the optical (g 475 -Z 85 o) and near- 
infrared (7iio-ffi6o) colors. The scatter about the linear calibration relations for the luminous red galaxies in 
the sample is approximately 0.10 mag, corresponding to a statistical error of 5% in distance. Our results imply 
that the distance to any suitably bright elliptical galaxy can be measured with this precision out to about 80 Mpc 
in a single-orbit observation with WFC3/IR, making this a remarkably powerful instrument for extragalactic 
distances. The calibration sample also includes much bluer and lower-luminosity galaxies than previously used 
for IR SBF studies, revealing interesting population differences that cause the calibration scatter to increase for 
dwarf galaxies. Comparisons with single-burst population models show that, as expected, the redder early-type 
galaxies contain old, metal-rich populations, while the bluer dwarf ellipticals contain a wider range of ages and 
lower metallicities than their more massive counterparts. Radial SBF gradients reveal that IR color gradients 
are largely an age effect; the bluer dwarfs typically have their youngest populations near their centers, while the 
redder giant ellipticals show only weak trends and in the opposite sense. Because of the population variations 
among bluer galaxies, distance measurements in the near-IR are best limited to red early-type galaxies. We 
conclude with some practical guidelines for using WFC3/IR to measure reliable SBF distances. 

Subject headings: distance scale — galaxies: clusters: individual (Virgo, Fornax) — galaxies: distances and 
redshifts — galaxies: elliptical and lenticular, cD — galaxies: stellar content 


1. INTRODUCTION 

1.1. Distance Measurements in the Big Picture 

Accurate distance measurement is central to both astro¬ 
physics and cosmology. Reliable distances are needed to con¬ 
vert observed properties of galaxies (fluxes and angular sizes) 

* Based on observations with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated with Program 
#11712, #11570, and #11691. 


into absolute physical quantities such as luminosities, masses, 
ages, star formation rates, and dynamical time scales. In the 
local universe where peculiar velocities are significant, the 
distance estimate is often a major source of uncertainty on 
these physical properties. For instance, a recent review ar¬ 
ticle on supermassive black holes (Kormendy & Ho 2013) 
notes that for many galaxies, the errors in the central black 
hole masses are dominated by the uncertainty in distance; yet, 
many authors neglect to include this important contribution to 
the uncertainty. 
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In the field of cosmology, the acceleration of the cosmic ex¬ 
pansion was first revealed by accurate distance measurements 
of Type la supernovae (Riess et al. 1998; Perlmutter et al. 
1999). The supernova distance estimates, combined with flat¬ 
ness constraints for the universe provided by the cosmic mi¬ 
crowave background (CMB) power spectrum, are primarily 
responsible for inaugurating a new era of “precision cosmol¬ 
ogy,” central to which is the conclusion that the mass-energy 
budget of the universe is dominated by “dark energy.” Now, 
with the exquisite constraints on the CMB power spectrum 
at z~1100 (e.g., Bennett et al. 2013; Planck Collaboration 
2015) and measurements of the baryon acoustic oscillation 
(BAO) scale at intermediate redshift (e.g., Eisenstein et al. 
2005; Blake et al. 2011), there has been a new effort to de¬ 
termine the local value of the Hubble constant Hq to a pre¬ 
cision of 1% (Riess et al. 2011; Freedman et al. 2012). As 
discussed in the foregoing references, this level of precision 
is required for firm simultaneous constraints on cosmic ge¬ 
ometry, the dark energy equation of state, and the number of 
neutrino species. 

Of course, the value of Hq has been controversial for 
decades, primarily owing to systematic calibration errors 
(e.g., Freedman & Madore 2010). Determining Hq with a to¬ 
tal uncertainty of no more than 1 % remains beyond the ability 
of any single distance measurement technique. In order to 
achieve the required level of precision, it is helpful to have 
multiple high-precision distance indicators to provide robust 
constraints on the contributions from systematic errors. 

The surface brightness fluctuation (SBF) method provides 
a measurement of the mean brightness of the red giant branch 
stars in an early-type galaxy even though individual stars can¬ 
not be resolved (Tonry & Schneider 1988). It was introduced 
as a way to estimate distances with ~ 10 % uncertainty out to 
about 20 Mpc from ground-based astronomical images (see 
Tonry et al. 1990). More recent applications have shown sub¬ 
stantial improvements in both the precision of the method and 
depth to which it can be applied (see the reviews by Blakeslee 
2012 and Fritz 2012), so that it has become one of a small 
number of methods capable of making a significant contribu¬ 
tion to the problem of constraining HqXo \%. In the follow¬ 
ing sections, we discuss these recent developments with the 
method and the need for a new calibration at near-IR wave¬ 
lengths. 

1.2. The Key Role of HST 

The latest revolution in our knowledge of the extragalactic 
distance scale (decreasing the uncertainty from nearly a factor 
of two to less than 10 %) has resulted primarily from obser¬ 
vations made with the Hubble Space Telescope (HST). Start¬ 
ing in the mid-1990s, HST was used to measure light curves 
for samples of Cepheid variable stars in late-type galaxies out 
to ~20 Mpc, mostly as part of the Key Project on the Dis¬ 
tance Scale (e.g., Freedman et al. 1994; Ferrarese et al. 1996; 
Kelson et al. 1996; Saha et al. 1996, 1997; Silbermann et al. 

1999) . The resulting Cepheid distance estimates were used 
to calibrate various secondary distance indicators and thereby 
derive the value of the Hubble constant Hq (Ferrarese et al. 
2000; Gibson et al. 2000; Sakai et al. 2000; Mould et al. 

2000) . The resulting value of Hq from the Key Project was 
72 ±8 kms“'Mpc“' (Freedman et al. 2001) where the total 
uncertainty includes both random and systematic contribu¬ 
tions. More recent Cepheid-based estimates are very similar 
to this value, but with reduced uncertainties of 3 to 4% (Riess 
et al. 2011; Freedman et al. 2012; Sorce et al. 2012). 


The excellent angular resolution and photometric stability 
of HST has also made it possible to measure SBF with far 
better precision and to much larger distances than was possi¬ 
ble from the ground. For instance, Jensen et al. (2001) cal¬ 
ibrated the SBF method for the near-IR F160W bandpass of 
the NICMOS NIC2 camera on HST (Thompson et al. 1999) 
and measured SBF distances to 16 galaxies beyond 40 Mpc, 
including the first SBF distances reaching beyond 100 Mpc. 
Stellar population effects on the NICMOS F160W SBF mag¬ 
nitudes were explored in detail by Jensen et al. (2003) using 
a larger sample of 65 galaxies. In general, the SBF absolute 
magnitude in a given bandpass depends on stellar population 
and must be calibrated using a population indicator, typically 
a broadband color. 

Following the installation of the Advanced Camera for Sur¬ 
veys (ACS) on HST, Mei et al. (2005b) produced the first cal¬ 
ibration of the SBF method for the ACS Wide Field Channel 
(ACSAVFC) using F475W (^ 475 ) and F850FP (zgjo) data for 
84 galaxies from the ACS Virgo Cluster Survey (ACSVCS; 
Cote et al. 2004). In that work, the zgso SBF measurements 
were calibrated for stellar population variations based on the 
observed (g 475 -Z 85 o) color; the resulting distances enabled the 
first clear resolution of the depth of the Virgo cluster and pro¬ 
vided constraints on its triaxial structure (Mei et al. 2007). In 
other studies based on ACSAVFC observations, Cantiello et 
al. (2005, 2007) measured multi-band SBF and color gradi¬ 
ents in 21 galaxies, and Biscardi et al. (2008) made the first 
optical SBF measurements beyond 100 Mpc. As part of the 
ACS Fornax Cluster Survey (ACSFCS; Jordan et al. 2007), 
Blakeslee et al. (2009) refined the SBF calibration for the 
ACS F850FP bandpass and determined the relative distance 
of the Virgo and Fornax clusters to a precision of 1.7%. Ad¬ 
ditional ACS/WFC SBF measurements and a new calibration 
for the F814W bandpass were published by Blakeslee et al. 
(2010). The launch of Wide Field Camera 3 with its powerful 
IR channel (WFC3/IR) has greatly increased the distance to 
which SBF measurements can be made within a single HST 
orbit. However, the method must first be calibrated for se¬ 
lected WFC3/IR passbands as done previously for NICMOS 
and ACS; this is the primary goal of the present work. 

1.3. SBF Measurements in the Infrared 

The development of new infrared detectors in the 1990s al¬ 
lowed researchers to successfully apply the SBF techniques 
to IR images for the first time (Fuppino & Tonry 1993; Pahre 
& Mould 1994; Jensen, Fuppino, & Tonry 1996). Because 
the SBF signal is dominated by the most luminous stars in a 
population, and these tend to be quite red for evolved galax¬ 
ies, SBF magnitudes are much brighter in the near-IR than at 
optical wavelengths. Additionally, extinction by dust (both in 
our Galaxy and in the target galaxy) is much lower at near- 
IR wavelengths. Depending on how it is distributed, dust can 
either reduce the fluctuation signal (as for a uniform screen 
of foreground dust in the Galaxy), or, more commonly, bias 
the fluctuation signal towards higher amplitudes and shorter 
SBF distances, as would occur if dusty regions were clumpy 
on scales comparable to the size of the point-spread function 
(PSF). Clumpy dust is often associated with recent star for¬ 
mation, and bright young stars seriously bias SBF measure¬ 
ments as well. The contrast between the fluctuations and other 
point-like sources (globular clusters and background galax¬ 
ies) is also higher in the near-IR bands, reducing yet another 
source of uncertainty in the SBF measurement. While the 
background in the IR is higher than at optical wavelengths. 
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the increased strength of the SBF signal more than compen¬ 
sates, especially from space, where the thermal background 
at 1.1 and 1.6 ^m is not significant. The benefit of the much 
lower background, combined with the excellent image qual¬ 
ity and a very stable PSF, usually makes near-IR SBF mea¬ 
surements with HST much more accurate than measurements 
from ground-based facilities. 

However, the calibration of the SBF magnitudes as a func¬ 
tion of stellar population is potentially more complicated in 
the near-IR. As noted above, the trend of SBF magnitude 
with galaxy color is used for calibrating the SBF distance 
measurements, both at optical and IR wavelengths. At op¬ 
tical wavelengths, the effects of age and metallicity variations 
on the SBF calibration relations are largely degenerate, but 
this degeneracy begins to break down in near-IR, and this 
can reveal interesting differences in the stellar populations 
of galaxies. Bluer elliptical and SO galaxies typically show 
signs of intermediate-age populations, and the asymptotic gi¬ 
ant branch (AGB) stars associated with those populations pro¬ 
duce brighter fluctuations (Jensen et al. 2003; Mieske, Hilker, 
& Infante 2003, 2006). The stellar population variations may 
therefore produce more scatter in IR SBF distance calibration, 
but the brightness of the fluctuations at these wavelengths 
make them measurable to much larger distances; thus, it is 
worth characterizing the behavior and limits of the calibration 
as well as possible. 

In this paper we report the results of a study to calibrate 
the IR SBF distance measurement technique using new SBF 
measurements in the FllOW (Tno) and F160W (//igo) band- 
passes of WFC3/IR on HST. The calibration sample includes 
16 galaxies spanning a wide range in galaxy luminosity and 
color. The following section describes the observations and 
sample in more detail. Section 3 discusses the data reductions 
and SBF measurements. New SBF calibrations in Tuo and 
are presented in Section 4, while implications of the SBF 
measurements for the galaxy stellar populations are discussed 
in Section 5. We provide our recommendations for measuring 
SBF distances with WFC3/IR in Section 6, before concluding 
with a summary. 

2. WFC3/IR OBSERVATIONS 

In order to calibrate the SBF method for WFC3/IR, we 
selected 16 early-type galaxies, eight in each of the Virgo 
and Fornax clusters, that already had high-quality ACS SBF 
measurements in zgso (g 475 -Z 85 o) colors (Blakeslee et al. 
2009). Table [T] lists the properties of the galaxies that were 
targeted in HST program GO-11712 (PI; J. Blakeslee). The 
galaxies were chosen to cover the full color range of the 
ACSVCS and ACSFCS samples so that the resulting cali¬ 
bration would be as generally applicable as possible. More¬ 
over, the NICMOS SBF calibration exhibited increased scat¬ 
ter for bluer galaxies (Jensen et al. 2003), as have ground- 
based /-band SBF measurements (Mieske et al. 2006); ex¬ 
ploring a broad color range should help us understand where 
the WFC3/IR calibration becomes less reliable. 

Each of these 16 galaxies were observed for one orbit, split 
approximately equally between the /no and //igo Alters, with 
four dithered exposures in each Alter and total exposure times 
varying with target visibility. The same dither pattern was 
used for each galaxy. 

Eor this study we also downloaded archival WEC3/IR //leo 
data forNGC 4258 (GO-11570, PI: A. Riess) and NGC 1316 
(GO-11691, PI: P. Goudfrooij). NGC 4258 is a late-type 
galaxy with a H 2 O masers in Keplerian orbits around a cen¬ 


tral black hole, enabling a geometric estimate of the distance 
(Greenhill et al. 1995; Miyoshi et al. 1995; Herrnstein et al. 
1999). While not ideal for SBE analysis, the importance of 
this galaxy to the absolute calibration of the extragalactic dis¬ 
tance scale makes it an important target worthy of a trial SBE 
measurement. NGC 1316 is an early-type SO galaxy in the 
Eornax cluster with extensive dust and signs of recent merg¬ 
ing. Although it is also a poor SBE candidate, it is a useful 
comparison galaxy for this study because it has hosted four 
type la supernovae. 

Tables |^and[^list the exposure times and sky brightnesses 
for the allthe Tuq and //igo observations used in this study. 
Additional information in these tables is discussed in the fol¬ 
lowing sections. 

3. SBE MEASUREMENTS 

The spatial fluctuations in the surface brightness of a 
smoothly distributed population of stars, as found in ellip¬ 
tical and lenticular galaxies, arise due to the Poisson statis¬ 
tics of the discrete stars making up the galaxy, even when 
the stars cannot be resolved or detected individually (Tonry 
& Schneider 1988). The SBE amplitude scales inversely with 
the square of the distance; nearby galaxies appear “bumpy” 
compared to more distant galaxies, where more stars are sam¬ 
pled by each resolution element and the \/N variation be¬ 
tween regions is therefore a smaller fraction of the total num¬ 
ber of stars. The fluctuations, which are dominated by the 
most luminous stars in a population, are blurTed by the PSE; 
additional contributions to the fluctuation signal arise from 
clumpy dust, globular clusters, background galaxies, and fore¬ 
ground stars. The process of making an SBE measurement 
consists of extracting and fitting the spatial Eourier power 
spectrum of the stellar fluctuations convolved with the PSE 
power spectrum and removing the contributions from extrane¬ 
ous sources. The resulting fluctuation power is used to com¬ 
pute the fluctuation magnitude. 

Procedures for measuring surface brightness fluctuations 
have been described in detail by several authors (e.g., Tonry 
et al. 1990, 1997; Blakeslee et al. 1997; Jensen et al. 1998; 
Mei et al. 2005a; Eritz 2012). The description here provides 
a concise overview of the process steps that are either unique 
to this study or are particularly relevant to the WEC3/IR SBE 
measurements. 

3.1. Data Reduction 

The first step in the SBE data reduction process involves 
producing a calibrated, combined, and background-subtracted 
image ready for further SBE analysis. We used the images 
from the HST archive reduced using the standard pipeline 
through the flat-fielding stage {fit files). Erom that point, we 
adopted a reduction procedure that differs from the standard 
pipeline. 

We combined the individual flat-fielded exposures using in¬ 
teger pixel shifts after fitting each image for background and 
identifying cosmic rays; in order to avoid introducing corre¬ 
lated noise between pixels, fractional pixel registration was 
not used. Using integer pixel shifts results in slightly lower 
spatial resolution in the combined image, but preserves the 
independence of noise from pixel to pixel, which is important 
for fitting the SBE power spectrum. 

In a second difference from the standard pipeline reduc¬ 
tion, the clean combined images were not corrected for the 
WEC3/IR geometrical distortion nor combined using astro- 
drizzle. Our analysis therefore includes the ~10% difference 
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TABLE 1 

Galaxy Properties 


Galaxy 

ClusteC 

Type'’ 

niB^ 

(mag) 

(arcsec) 

Ms” 

(mag) 

Rj 

(kpc) 

Ay8 

(mag) 

Alt id'’ 

IC 1919 

F 

dSO 

13.5 

21.2 

-18.1 

2.06 

0.010 

FCC 43 

IC 2006 

F 

El 

12.2 

18.9 

-19.4 

1.84 

0.008 


NGC 1344 

F 

E5 

11.3 

33.2 

-20.3 

3.22 

0.013 

NGC 1340 

NGC 1374 

F 

EO 

11.9 

29.2 

-19.7 

2.83 

0.010 

ECC 147 

NGC 1375 

F 

SO 

13.6 

13.2 

-18.0 

1.28 

0.010 

ECC 148 

NGC 1380 

F 

SO 

11.3 

30.3 

-20.3 

2.94 

0.012 

ECC 167 

NGC 1399 

F 

EO 

10.6 

114.6 

-21.0 

11.11 

0.009 

ECC 213 

NGC 1404 

F 

E2 

10.9 

22.9 

-20.7 

2.22 

0.008 

ECC 219 

IC 3025 

V 

dSO 

14.8 

9.8 

-16.4 

0.78 

0.015 

VCC21 

IC 3032 

V 

dE2 

14.7 

9.1 

-16.6 

0.73 

0.026 

VCC33 

IC 3487 

V 

dE6 

14.8 

9.3 

-16.4 

0.74 

0.015 

VCC 1488 

IC 3586 

V 

dSO 

14.5 

20.2 

-16.7 

1.61 

0.032 

VCC 1695 

NGC 4458 

V 

El 

12.9 

20.4 

-18.2 

1.63 

0.017 

VCC 1146 

NGC 4472 

V 

E2 

9.3 

210.8 

-21.9 

16.87 

0.016 

VCC 1226 

NGC 4489 

V 

SO 

12.8 

37.1 

-18.4 

2.97 

0.020 

VCC 1321 

NGC 4649 

V 

E2 

9.8 

98.1 

-21.4 

7.85 

0.019 

VCC 1978 


^ Cluster: V for Virgo and F for Fornax. 

^ Morphological type from the ACSVCS (Cote et al. 2004) and ACSFCS (Jordan et al. 2007). 
Apparent B-band magnitude (Vega). 

^ Effective radius in arcseconds, determined from the ACS and/or SDSS imaging (Ferrarese et al. 
2006; Chen et al. 2010; P. Cote, priv. comm.). 

^ Absolute B magnitude (Vega), con'ected for Galactic extinction, and assuming distances of 16.5 
and 20.0 Mpc for galaxies in Virgo and Fornax, respectively. 

^ Effective radius in kpc, assuming same Virgo and Fornax distances as above. 

§ Galactic extinction (Vega mag) in 7-band, from Schlafiy & Finkbeiner (2011). 

^ Alternate names from Virgo and Fornax Cluster Catalogs (Binggeli et al. 1985; Ferguson 1989), 
or alternate NGC designation in the case of NGC 1344. 


in plate scale between the x and y axes, causing our images 
to appear somewhat narrower horizontally than they do on the 
sky (Fig. 0. The SBF procedure involves taking the spatial 
Fourier power spectrum, and correlated noise between pix¬ 
els resulting from fractional pixel shifts and interpolated pixel 
values when correcting for geometrical distortion can produce 
a slope in the white noise component of the power spectrum. 

Finally, we chose not to apply the correction to pixel size 
in the y-axis of WFC3/IR images. The pixel map correction 
usually used for WFC3/1R corrects for PSF variations but is 
inappropriate for extended objects. The WFC3/1R focal plane 
is tilted, and the size (area) of the pixels on the sky varies 
by ~8% from the center to the upper and lower edges of 
the frame (Kalirai et al. 2010). When the data are divided 
by the flat held image (as is done for the fit hies in the HST 
archive), the varying pixel sensitivity is removed, effectively 
making the pixels equally sensitive to uniform illumination. 
Flattening the images in this way creates a variation in sensi¬ 
tivity for point sources from the center to the top and bottom 
edges that is usually corrected in the pipeline data reduction 
process using astrodrizzle. Since we are interested in accu¬ 
rately measuring the surface brightnesses of the galaxies and 
avoiding correlated noise between pixels, we chose not to cor¬ 
rect for pixel size variation in our SBF data reduction process. 
The practical effect of this decision is that the PSF at the ex¬ 
treme upper and lower regions of the image does not match 
the PSF near the center. We have carefully chosen PSF refer¬ 
ence stars from the same vertical region of the field of view 
as the galaxies being analyzed (usually very near the center) 
to avoid any systematic offset between the PSF photometric 
normalization or power spectrum shape and the galaxy fluc¬ 
tuations. The chosen PSF stars were all unresolved, isolated 
from other bright objects, and much brighter than the globular 
cluster population. 


3.2. Measuring Fluctuations 

SBF measurements are made by htting the Fourier spatial 
power spectrum of the stellar fluctuations in a given region 
of a galaxy with the normalized power spectrum of the PSF. 
There are several steps to produce the spatial power spec¬ 
trum: (i) the background level is estimated and subtracted; 

(ii) extraneous objects (globular clusters, background galax¬ 
ies, and dusty regions) are identified, measured, and masked; 

(iii) a smooth isophotal model is fitted to the galaxy profile 
and subtracted; (iv) isolated bright stars are extracted and used 
to determine the power spectrum of the PSF; (v) the Fourier 
power spectra are computed for the stellar fluctuations, the 
masked galaxy profile, and the PSF; and (vi) the power spec¬ 
trum of the data is fitted and normalized to determine the SBF 
power in units of flux, with the power from undetected globu¬ 
lar clusters and galaxies subtracted, from which a fluctuation 
magnitude is computed in the established way (e.g., Tonry & 
Schneider 1988; Jensen et al. 1998; Mei et al. 2005a). 

It is important to accurately measure and subtract the IR 
background before computing fluctuation magnitudes. Esti¬ 
mating the background level was done in an iterative process 
of cross-checking values measured in different ways. For 
the small galaxies, we measured the flux in the comers of 
the frames. For the larger ellipticals, we also determined the 
background using the best fit to a r'/'* profile for each galaxy. 
These estimates were then compared to a measurement of the 
background made by iteratively computing a smooth model 
for each galaxy and looking at the residual background in the 
held of view. By adjusting the sky level offset, we optimized 
the galaxy models such that the residual background was not 
systematically positive or negative. The difference between 
the sky values determined using the different methods was 
used as an estimate of the uncertainty in the sky level, and 
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TABLE 2 

7iio SBF Measurements 


Galaxy 

Exposure 

(sec) 

Background 
(AB mag/arcsec^) 

Annulus 

(arcsec) 

{gal/sky)‘‘ 

average 

gal/sky 

range 

SBFS/A?'’ 

(Po-Pr)/Pi 

mno" 

(AB mag) 

Fornax 

IC 1919 

997 

21.74 

8-33 

1.3 

0.9-2.6 

33 

28.31 ±0.05 

IC 2006 

997 

22.40 

4-33 

12 

6-58 

43 

28.63 ±0.05 

NGC 1344 

997 

22.10 

4-33 

19 

11-85 

83 

28.37 ±0.05 

NGC 1374 

997 

22.02 

4-33 

10 

5-48 

47 

28.55 ±0.06 

NGC 1375 

1197 

22.64 

4-33 

5.8 

3.1-27 

59 

28.19 ±0.05 

NGC 1380 

1197 

22.22 

6-33 

29 

19-124 

49 

28.51 ±0.04 

NGC 1399 

1197 

22.15 

4-33 

41 

23-194 

52 

28.79 ±0.05 

NGC 1404 

997 

22.04 

4-33 

32 

17-173 

47 

28.68 ±0.06 

Virgo 

IC 3025 

997 

21.84 

8-17 

0.8 

0.8 

17 

28.36 ±0.06 

IC 3032 

997 

21.81 

8-17 

0.8 

0.8 

34 

27.95 ±0.07 

IC 3487 

997 

21.84 

4-33 

0.4 

0.2-2.5 

27 

28.16±0.10 

IC 3586 

997 

22.17 

4-33 

0.8 

0.5-3.6 

34 

28.05 ±0.08 

NGC 4458 

997 

22.12 

4-67 

4.4 

0.5-22 

47 

28.01 ±0.05 

NCG 4472 

748 

22.10 

4-67 

68 

17-249 

40 

28.38 ±0.06 

NGC 4489 

997 

22.23 

8-67 

4.5 

0.6-8 

58 

27.81 ±0.06 

NGC 4649 

997 

22.10 

4-67 

61 

14-250 

43 

28.44 ±0.06 

Supernova host 

NGC 1316 

1396 

22.08 

33-67 

15 

8-34 

166 

28.26 ±0.07 


^ Weighted average ratio of the galaxy surface brightness to sky background surface brightness within the measurement 
region. 

^ Weighted average ratio of the SBF fluctuation power to the white noise component Pi of the spatial power spectrum. 
SBF magnitudes have been corrected for Galactic extinction using values from Schlafly & Finkbeiner (2011). 


TABLE 3 

Hieo SBF Measurements 


Galaxy 

Exposure 

(sec) 

Background 
(AB mag/arcsec^) 

Annulus 

(arcsec) 

{gal/sky)‘‘ 

average 

gal/sky 

range 

SBFS/A?'’ 

{Po-Pr)lPi 

mi60^ 

(AB mag) 

Fornax 

IC 1919 

997 

21.90 

8-33 

1.8 

1.3-3.8 

43 

27.55 ±0.06 

IC 2006 

997 

22.13 

4-33 

12 

6-59 

59 

27.86 ±0.04 

NGC 1344 

997 

21.90 

4-33 

20 

12-91 

77 

27.52 ±0.05 

NGC 1374 

748 

21.95 

4-33 

12 

6-59 

45 

27.88 ±0.05 

NGC 1375 

1197 

22.36 

4-33 

5.5 

2.9-27 

62 

27.42 ±0.05 

NGC 1380 

1197 

21.95 

6-67 

29 

19-120 

71 

27.80 ±0.05 

NGC 1399 

1197 

21.95 

4-33 

46 

25-220 

52 

28.03 ±0.04 

NGC 1404 

997 

22.07 

4-33 

43 

22-235 

50 

27.94 ±0.06 

Virgo 

IC 3025 

997 

21.76 

8-17 

0.3 

0.3 

27 

27.72 ±0.07 

IC 3032 

997 

21.86 

8-17 

1.0 

1.0 

17 

26.94 ±0.09 

IC 3487 

997 

21.75 

4-33 

0.4 

0.2-2.7 

27 

27.42 ±0.09 

IC 3586 

997 

22.01 

4-33 

0.9 

0.5-3.8 

38 

27.26 ±0.04 

NGC 4458 

997 

21.94 

4-67 

4.7 

0.5-23 

56 

27.35 ±0.05 

NCG 4472 

748 

21.83 

4-67 

71 

17-258 

47 

27.59 ±0.05 

NGC 4489 

997 

22.07 

8-67 

4.8 

0.6-8 

67 

27.10±0.05 

NGC 4649 

997 

21.90 

4-67 

68 

15-283 

67 

27.63 ±0.05 

Supernova host 

NGC 1316 

2796 

21.85 

33-67 

21 

9-36 

111 

27.35 ±0.07 

Maser host 

NGC 4258 

2012 

21.91 

irreg. 

10 

47 

47 

25.30 ±0.06 


^ Weighted average ratio of the galaxy surface brightness to sky background surface brightness within the measurement 
region. 

^ Weighted average ratio of the SBF fluctuation power to the white noise component of the spatial power spectrum. 

SBF magnitudes have been corrected for Galactic extinction using values from Schlafly & Finkbeiner (2011). 
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IC 1919 1.163 


f 


NGC 1375 1.256 


IC 3025 0.919 


IC 3032 1.006 



“ NGC 1344 1.319 NCG 1380 1.391 NGC 1404 1.471 




NGC 1374 1.375 IC 2006 1.409 NGC 1399 1.490 



IC 3487 1.068 NGC 4458 1.236 NGC 4472 1.514 



IC 3586 1.118 NGC 4489 1.257 NGC 4649 1.554 


Fig. 1.— Fornax (top two rows) and Virgo cluster galaxies (bottom two rows). Each image is displayed with background subtracted and using the same upper 
and lower limits and logarithmic stretch, within the same field of view 100 arcsec on a side. The bluer galaxies are on the left and the redder giant ellipticals are 
on the right. Measured values of (^ 475 —^ 850 ) (AB) are shown next to the galaxy labels. 


Infrared SBF Distances Using WFC3 


7 


|iiii|iiii|iiii|iiii|it 

IC1919 ■ 




iii|iiii|iiii|iiii|i(- 

N1404 ^ 


a> 

o 

04 

(30 

o 



3.0 

2.0 

1.0 

3.0 

2.0 

1.0 

3.0 

2.0 

1.0 

3.0 

2.0 — 

1.0 

liiiiliiiiliiiiliiiilirHiiiiliiiiliiiiliiiilirniiiiliiiiliiiiliiiilirHiiiiliiiiliiiiliiiili 
0 50 100 150 200 0 50 100 150 200 0 50 100 150 200 0 50 100 150 200 

Wavenumber (fc) 

Fig. 2.— Fits to the spatial power spectra for the 16 galaxies in Fornax and Virgo (superimposed white line). The dashed lines indicate the scaled PSF power 
spectrum and the white noise (flat) components. 


the SBF analysis was repeated to determine the uncertainty in 
fluctuation magnitude due to sky level uncertainty. The SBF 
magnitude is normalized by the mean galaxy brightness at the 
location of the SBF measurement. For the smaller galaxies 
with the lowest surface brightness (see Fig.[T]i, the background 
level was most accurately measured. For the largest galaxies 
that extend beyond the limits of the field of view, the back¬ 
ground was most difficult to measure accurately, but its influ¬ 
ence on the SBF measurement was also minimal. The uncer¬ 
tainty in the SBF magnitude is therefore relatively insensitive 
to the uncertainty in the background measurement. 

Typical background levels at /no were 1.3 e“s“*pix“', or 
22.10±0.23 AB magarcsec"^; at i/ieo we measured 0.67 
e“s“*pix“\ or21.95±0.15 AB magarcsec"^ (see Tables|^and 
[^. These sky values were found to be consistent with pub¬ 
lished empirical measurements of the HST IR background, 
which is dominated by scattered zodiacal light in the Solar 
System (Pirzkal 2014). The /no background levels were also 
affected by the diffuse upper-atmosphere He emission line 
at 1.083 /im (Brammer et al. 2014). While very few indi¬ 
vidual exposures were strongly affected, the residual back¬ 
ground variation among the data is larger at 7no than at //leo. 
The background variability between observations of different 
galaxies (see Tables and is not directly related to the 
uncertainty in the baSground measurement for a particular 
galaxy, which was determined independently for each galaxy. 

After the background had been subtracted, we then iden¬ 


tified and masked any non-fluctuation point sources in the 
held of view. Background galaxies and globular clusters were 
identified and their brightnesses measured using the SExtrac- 
tor software (Bertin & Arnouts 1996) and adopting aperture 
corrections taken from the WFC3 instrument webpagej^ We 
started by making an initial fit to the galaxy and subtracting 
it to make it easier for SExtractor to identify and measure all 
the compact objects in the held (we used a noise model for 
SExtractor that accounts for the subtracted galaxy and pre¬ 
vents the software from confusing the stellar fluctuations with 
globular clusters). Using the SExtractor output, we then cre¬ 
ated a fit to the luminosity functions of the globular clusters 
and background galaxies. Eor most of the galaxies, we used a 
Gaussian of width 1.2 mag to fit the luminosity functions; we 
found that a width of 1.35 mag fit better for NGC 1399 and 
NGC 4472. We adopted a luminosity function peak absolute 
magnitude of My=-1 AQ, My=-8.26, and Mh=- 8.29 (e.g., 
Harris 2001; Erogel et al. 1978). The background galaxies are 
assumed to follow a power law distribution with power law 
slope of 0.25, which for most of our observations results in a 
normalization of about one to two galaxies per square arc sec 
at 34.50 mag AB (Retzlaff et al. 2010; Windhorst et al. 2011). 
These fits allowed us to integrate the contribution to the SBE 
signal fainter than the completeness limit and correct the final 
SBE magnitude accordingly. A mask was then created from 
the list of objects that removed all objects brighter than the 

^ http://www.stsci.edu/hst/wfc3/phot_zp_lbn 
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limiting cutoff magnitude. 

The next step was to fit a final smooth model to the galaxy 
profile with the sky and point sources removed. We used an 
iterative procedure of fitting elliptical isophotes to the galaxy 
surface brightness profile, allowing the procedure to adjust 
the centers, ellipticities, and orientations of the elliptical aper¬ 
tures. This model galaxy was later used to normalize the fluc¬ 
tuation signal. 

With the sky background and mean galaxy surface bright¬ 
ness removed and contaminating point sources masked, we 
then computed the Fourier transform and spatial power spec¬ 
trum of the data in circular annuli (see Tables and for 
the sizes of the annular regions). The measurements were re¬ 
peated using elliptical annuli for the six most elongated galax¬ 
ies with apparent radial gradients in their SBF magnitudes. 
The power spectrum was normalized by the mean galaxy lu¬ 
minosity so that the fluctuation amplitude should be the same 
in each annulus. The purpose of measuring fluctuations in 
multiple annuli was to look for consistency between regions 
with varying surface brightness, globular cluster population, 
and distance from the center. It also allowed us to measure 
the radial gradients in fluctuation ampl itude, and thus stellar 
population, as described in Section [ 53 ] 

The measured fluctuation power spectrum is a convolution 
of the pixel-to-pixel variation in the number of stars and the 
PSF of the instrument. We therefore require a robust measure¬ 
ment of the PSF power spectrum to determine the fluctuation 
amplitude. We extracted isolated bright stars from the central 
region of the detector field of view and computed their power 
spectra. Because PSF stars are not all uniformly centered on 
the pixels, there is some variation that naturally arises in the 
PSF power spectra. To determine which PSF stars best fit¬ 
ted the power spectrum for a particular galaxy, we combined 
the cleanest and brightest PSF stars from several observations, 
and then repeated the SBF measurements using a variety of 
PSF stars. We computed the uncertainty in SBF magnitude 
attributable to the PSF variations by determining the range of 
plausible PSF fits based on the shape of the power spectra 
and the quality of the fits. We then used a common compos¬ 
ite PSF to measure SBF consistently in all the galaxies (the 
same dither pattern was used for all the observations). As an 
added check on the PSF uniformity and fit quality, we also 
fitted the observed power spectra to the “Tiny Tim” PSF mod- 
el^(Krist et al. 2011). The Tiny Tim model PSF for each 
filter was convolved with Gaussians of various widths to con¬ 
struct a library of model PSFs to provide better matches to the 
data, which had been combined using integer p ixel offsets and 
no geometrical distortion corrections (see Sec. |3.1| l. The SBF 
magnitudes computed using the Tiny Tim model PSFs were 
then compared to those derived from the combined empirical 
PSFs and the range of plausible fits was used to determine the 
uncertainty due to PSF variations (typically 0.04 mag). 

Fluctuation magnitudes were computed by fitting the nor¬ 
malized PSF power spectrum to the galaxy power spectrum 
P{k) = P()E(k)+Pi, where E(k) is the expectation power spec¬ 
trum, which includes the smooth galaxy profile and the com¬ 
bined annular region and external object mask, all convolved 
with the normalized PSF power spectrum (see Fig. 13 for 
power spectra and fit components). The fits excluded the 
lowest wavenumbers k < 10, which are affected by large-scale 
galaxy and sky subtraction errors. Pi is the white noise com¬ 
ponent, which is flat for uncorrelated pixel-to-pixel noise. The 

^ http://tinytim.stsci.edu/cgi-bin/tinytimweb.cgi 


scale factor used to best match the data corresponds to the flux 
in SBF power in units of e“s“'. It is then straight-forward to 
compute the fluctuation magnitude: 

m = -2.5 log(Po - Pr)+Ml 

where Pq is the fluctuation power and P^ is the contribution 
from point sources fainter than the completeness limit. Mi 
is the zero point for the filter-tdetector (26.8223 AB for Jhq 
and 25.9463 AB for //leo)- The AB magnitude is 0.7595 mag 
larger than the Vega equivalent at /no, and 1.2514 mag larger 
than the Vega magnitude at //leo- AB magnitudes and colors 
were corrected for Galactic extinction using the values pub¬ 
lished by Schlafly & Finkbeiner (2011). The (g 475 -Z 85 o) col¬ 
ors from Blakeslee et al. (2009) were adjusted to make them 
consistent with the Schlafly & Finkbeiner extinction values. 
Individual sky background levels, e^osure times, and SBF 
S/N values are listed in Tables [^ant^ 

Absolute fluctuation magnitudes Muq and Migg (Table l^i 
were then computed using both individual distance modulus 
measurements and average cluster distances of (m-M) =31.51 
mag for the Fornax cluster and 31.09 for Virgo (all from the 
optical ACS SBF measurements of Blakeslee et al. 2009). The 
radial extent of each cluster (0.053 mag for Fornax and 0.085 
mag for Virgo) was adopted as the uncertainty on the aver¬ 
age cluster distance moduli. The colors originally reported 
by Blakeslee et al. (2009) were recomputed for the apertures 
used in this study (circular and elliptical). 

There are several sources of uncertainty in the SBF mea¬ 
surement procedure that we quantified by exploring the range 
of input parameters, as we did for the uncertainty due to the 
PSF fit. Average values (and ranges) of the uncertainties we 
measured and adopted for the final SBF measurements are 
listed in Table [3 Not all sources of uncertainty are com¬ 
pletely independent—^residual errors in sky subtraction can 
affect the galaxy or PSF fit, for example—so the total uncer¬ 
tainties listed in Table are not a simple quadrature addition 
of all sources listed in Table |3 we estimated the fraction of 
the power spectrum fit (Pq) uncertainty that results from the 
PSF fit and sky subtraction separately before adding all inde¬ 
pendent sources of error in quadrature. The distance modulus 
un(^rtainties were includedjri the individual-distance values 
of M. The cluster-distance M values include the cluster dis¬ 
tance dispersion values from Blakeslee et al. (2009) in the 
total uncertainty. 

4. ANALYSIS 

4.1. Calibration of the WFC3/IR SBF Distance Scale 

The SBF signal can be detected with HST in modest 7no 
and //i6o exposures (an orbit or less) out to ^100 Mpc (e.g., 
Jensen et al. 2001). To take full advantage of WFC3/IR ob¬ 
servations of early-type galaxies collected for a variety of 
purposes and measure accurate distances, we calibrated the 
WFC3/IR 7iio and //igo SBF distances by fitting the abso¬ 
lute fluctuation magnitudes M as a function of both optical 
(g 475 -Z 85 o) and IR (7iio-77i6o) colors. Determining the M val¬ 
ues requires us to adopt a distance modulus {m-M) for each 
galaxy. We used the Z 850 SBF distance moduli measured using 
ACS (Blakeslee et al. 2009). These measurements provide a 
consistent distance reference accurate to < 0.08 mag for most 
of the galaxies, although the calibration becomes systemati¬ 
cally less certain for (^ 475 -^ 850 ) < 1 -05 mag. We also adopted 
(^ 475 “Z 85 o) values from the same ACS data, supplementing 
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TABLE 4 

Absolute SBF Magnitudes 


Galaxy 

(m-MY 

(g475“Z850)’’ 

UllO“^^16o) 



^160“ 


Fornax 

IC 1919 

31.485±0.073 

1.163 ±0.037 

0.223 ±0.007 

-3.18 ±0.09 

-3.20 ±0.08 

-3.94 ±0.09 

-3.96 ±0.08 

IC 2006 

31.525 ±0.086 

1.409±0.013 

0.263 ±0.013 

-2.90 ±0.10 

-2.88 ±0.07 

-3.67 ±0.10 

-3.65 ±0.07 

NGC 1344 

31.603 ±0.068 

1.319±0.007 

0.260 ±0.008 

-3.23 ±0.08 

-3.14±0.07 

-4.08 ±0.08 

-3.99 ±0.07 

NGC 1374 

31.458 ±0.070 

1.375±0.011 

0.251 ±0.014 

-2.91 ±0.09 

-2.96 ±0.08 

-3.58 ±0.09 

-3.63 ±0.07 

NGC 1375 

31.500 ±0.072 

1.256 ±0.029 

0.224 ±0.016 

-3.31 ±0.09 

-3.32 ±0.07 

-4.08 ±0.09 

-4.09 ±0.07 

NGC 1380 

31.632±0.075 

1.391 ±0.007 

0.286 ±0.006 

-3.12 ±0.09 

-3.00 ±0.07 

-3.83 ±0.09 

-3.71 ±0.07 

NGC 1399 

31.596±0.091 

1.490 ±0.005 

0.302 ±0.012 

-2.81±0.10 

-2.72 ±0.07 

-3.57±0.10 

-3.48 ±0.07 

NGC 1404 

31.526 ±0.072 

1.471 ±0.006 

0.292 ±0.005 

-2.85 ±0.09 

-2.83 ±0.08 

-3.59±0.10 

-3.57 ±0.08 

Virgo 

IC 3025 

31.421 ±0.130 

0.919 ±0.074 

0.183±0.019 

-3.06 ±0.14 

-2.73 ±0.11 

-3.70±0.15 

-3.37±0.11 

IC 3032 

30.886 ±0.133 

1.006 ±0.030 

0.184±0.036 

-2.94 ±0.15 

-3.14±0.11 

-3.95±0.16 

-4.15±0.13 

IC 3487 

31.053 ±0.134 

1.068 ±0.060 

0.132±0.023 

-2.89±0.17 

-2.93 ±0.13 

-3.63±0.16 

-3.67±0.13 

IC 3586 

31.093 ±0.080 

1.118±0.040 

0.188 ±0.007 

-3.04 ±0.11 

-3.04 ±0.11 

-3.83 ±0.09 

-3.83±0.13 

NGC 4458 

31.063 ±0.070 

1.236 ±0.049 

0.237 ±0.026 

-3.05 ±0.08 

-3.08 ±0.10 

-3.71 ±0.09 

-3.74 ±0.10 

NCG 4472 

31.116±0.075 

1.514±0.006 

0.291 ±0.013 

-2.74 ±0.10 

-2.71 ±0.10 

-3.53 ±0.09 

-3.50±0.10 

NGC 4489 

30.935 ±0.069 

1.257±0.014 

0.226 ±0.013 

-3.13±0.09 

-3.28±0.10 

-3.84 ±0.08 

-3.99±0.10 

NGC 4649 

31.082 ±0.079 

1.554 ±0.006 

0.311 ±0.013 

-2.64 ±0.10 

-2.65 ±0.11 

-3.45 ±0.09 

-3.46 ±0.10 

Supernova host 

NGC 1316 

31.606 ±0.065 

1.374 ±0.007 

0.272 ±0.007 

-3.35±0.10 

-3.25 ±0.09 

-4.26 ±0.10 

-4.16 ±0.09 

Maser host 

NGC 4258 

29.404 ± 0.048 

1.361 ±0.044 




-4.10 ±0.08 



Note. — All magnitudes are on the AB system and extinction corrected. 

“ Blakeslee et al. (2009) except NGC 4258, Humphreys et al. (2013). 

** Galaxy colors from Blakeslee et al. (2009) have been updated to match the apertures used in this study, and corrected for extinction 
using Schlafly & Finkbeiner (2011). 

M computed using the individual distance moduli shown in the second column. 

^ M computed using average cluster distances of 31.51 mag for Fornax and 31.09 for Virgo. Uncertainties on M include the cluster depths 
of 0.053 mag (Fornax) and 0.085 mag (Virgo). 


TABLE 5 

Average Uncertainties 


Source 

cr (mag) 

Range (mag) 

Power spectrum fit Pq 

0.06 

0.04-0.10 

PSF fit 

0.04 

0.026-0.067 

Background subtraction 

0.01 

0.002-0.062 

Galaxy subtraction 

0.007 

0 .001-0.022 

(g475-Z85o) color uncertainty 

0.025 

0.005-0.074 

color uncertainty 

0.013 

0.002-0.036 

Distance modulus (individual) 

0.086 

0.068-0.134 

Distance modulus (Virgo) 

0.085 

0.085 

Distance modulus (Fornax) 

0.053 

0.053 

Total statistical M uncertainty 

0.10 

0.08-0.17 


the published values from Blakeslee et al. (2009) with up¬ 
dated color measurements made using the original images in 
annuli that matched our IR observations. The ACS colors and 
distance moduli used to calibrate the IR SBF distance scale 
are listed in Table Extinction-corrected (7iio-^fi6o) colors 
were measured using our WFC3/IR images. 

As in optical bandpasses, the intrinsic luminosity of IR fluc¬ 
tuations varies with galaxy color. The SBF amplitude is sen¬ 
sitive to the brightness of the most luminous stars in a pop¬ 
ulation, and bluer galaxies with a significant component of 
young or intermediate-age stellar populations have luminous 
AGB stars that enhance the SBF signal (e.g., Jensen et al. 
2003). An accurate IR SBF distance calibration must take 
into account the brightening of fluctuations at intermediate 
and bluer colors. As in the previous F160W NICMOS cali¬ 
bration (Jensen et al. 2003), a linear fit to the red end of the 
sample was found to best represent the SBF calibration for 
distance measurements of typical giant elliptical galaxies. To 
determine the best slope of M as a function of galaxy color 


(see Figs. and 1^, we adopted an iterative procedure that 
takes into account the uncertainties in both the M and color 
axes. We started by making an initial approximate fit ignoring 
the color uncertainties (i.e., using a standard least-squares ap¬ 
proach). We computed the eiTor ellipse for each point and the 
distance from the fitted line in units of the combined x and y 
uncertainties. Because the zgso SBF distances are also a func¬ 
tion of (g475-Z85o)^ there is a small coiTelation between the x 
and y-axis uncertainties when fitting the individual-distance 
M values vs. (g 475 -Z 85 o)- We included the rotation of the er¬ 
ror ellipse for that subset of the calibration fits. We then iter¬ 
atively adjusted the fit coefficients to minimize the combined 
difference between the line and the data points, and then com¬ 
puted the rms in M. This procedure was repeated for each of 
the filters for both (g 475 -Z 85 o) and (/no-J^ieo)^ and for the two 
sets of M values derived from individual and common clus¬ 
ter distances. The coefficients and rms scatter for each fit are 
listed in Tables |6] and I?] 

The intrinsic scatter in the optical SBF distance measure¬ 
ments is comparable to the Virgo and Fornax cluster depths 
as estimated by Blakeslee et al. (2009) (0.053 mag for For¬ 
nax, 0.085 mag for Virgo). Depending upon the location of 
the galaxy within the cluster, the error in the estimated dis¬ 
tance may be larger when adopting the cluster mean or when 
using the individual galaxy distance. If the typical distance 
measurement error (including intrinsic scatter about the stel¬ 
lar population calibration) is less than the magnitude of the 
scatter from cluster depth, then it makes sense to calibrate the 
IR measurements using individual optical SBF distances. If 
the optical SBF distance measurement errors dominate, then 
it would be better to average all the optical SBF measure¬ 
ments and use a common cluster distance to calibrate the IR 
measurements. In the case of Fornax, the individual distance 
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Fig. 3.— Fits to the absolute fluctuation magnitudes M as a function of (^ 475 —Z 85 o) color. We also plot M for NGC 4258 for comparison. The calibration shown 
on the left is based on distances derived using the individual z-band SBF distances. The panel on the right shows the calibration using average cluster distances 
for Virgo and Fornax. Red points indicate //160 measurements (top set) and blue symbols are the ^uo measurements (lower set). The dashed calibration lines 
show the quadratic fit and the linear calibration including IC 1919; the color of IC 1919 is intermediate between the other blue dwarf galaxies and the elliptical 
galaxies in our sample. 


errors are larger than the scatter from cluster depth; the situ¬ 
ation is less clear in Virgo, especially for the bluer galaxies. 
Table [^presents calibrations computed using both individual 
distance moduli and average cluster distances of 31.09 mag 
and 31.51 mag for Virgo and Fornax, respectively (Blakeslee 
et al. 2009). The linear calibrations using cluster distances are 
not significantly different, particularly for the bright ellipticals 
located near the centers of the clusters, but the rms scatter is 
somewhat lower using the cluster distances. 

The downturn in absolute SBF magnitude at the blue end 
is dominated by the four dwarf galaxies in Virgo. The four 
bluest galaxies have the largest color uncertainties and ranges, 
the lowest SBF S/N ratios, the lowest galaxy brightness com¬ 
pared to the sky, and the largest radial gradients in fluctuation 
amplitude. They also show evidence of a wide range in stellar 
population age and metallicity, as manifested in their large ra¬ 
dial color and fluctuation magnitude gradients (the error bars 
shown in Figures]^ andare larger than the measurement er¬ 
rors, and include me range o f val ues due to radial gradients 
as described below in Section [53] l. These four Virgo galaxies 
were excluded from the linear calibration fits; IC 1919 was 
excluded from the (g 475 -Z 85 o) fits as well. Because the fluc¬ 
tuation amplitude is significantly lower in the bluest galaxies 
in this sample, we also present an alternative second order 
polynomial distance calibration fit that can be used for lower 
accuracy distance measurements of bluer galaxies (Figs.|^and 
1^. The quadratic fits were computed using the same iterative 
procedure that was used to make the linear fits, as described 
above. Higher-order fits are not justified given the sample 
size, measurement uncertainties, and large population varia¬ 
tions between blue dwarf elliptical galaxies. The linear cali¬ 
bration is not useful for the bluer galaxies; they have too much 
population variation for SBF to be generally useful as a dis¬ 
tance indicator. The quadratic calibration may be used to get 
approximate distances for bluer galaxies when necessary. 

Calibration coefficients are shown in Tables |6] and |7] The 


coefficients are defined as follows: 

M = a + b[{g-z)-lA]+c[{g-z)-lAf (1) 

M = a +b[(J-H)-0.27] +c[(J-H)-0.27 f (2) 

where c=0 for the linear fits. 

For (g475-Z85o)>l-2, researchers measuring distances 
should use the appropriate linear calibrations centered at the 
mean galaxy color as follows: 

Miio = (-2.946±0.015)-i-(2.16±0.15)[(g-z)-1.4] (3) 
Mm = (-3.699±0.028)-r(2.13±0.27)[(g-z)-1.4] . (4) 

The quadratic fits may be used for bluer galaxies, bearing 
in mind that the intrinsic scatter between blue galaxies is 
large. Quadratic fits are not shown for M computed using 
cluster distances because two of the bluest Virgo galaxies 
have individual optical SBF distance moduli that differ from 
the mean cluster modulus for Virgo by -0.20 and h- 0.33 mag 
(IC 3032 and IC 3025, respectively). The differences are 2.4 
and 3.9 times larger than the Virgo cluster depth of 0.085 mag 
(Blakeslee et al. 2009); these galaxies are probably outside the 
Virgo cluster core and should not be included in a calibration 
based on mean cluster distances. 

If (^ 475 “Z 85 o) colors are not available, distances may be 
computed using the {J\m-Hm) color instead. The scat¬ 
ter in the calibration with (JuQ-Hm) is larger than with 
(g 475 -Z 85 o) because the color range spanned is much smaller. 
For (7 iio-//i6o)> 0.2, the following relations centered at the 
mean galaxy color should be used: 

Mno = (-2.964±0.032)-i-(6.7±0.9)[(7-//)-0.27] (5) 
Mi6o = (-3.718±0.035)-i-(7.1± 1.1)[(7-//)-0.27] . (6) 

For magnitudes on the Vega system, subtract 0.7595 mag from 
7iio AB and 1.2514 mag from Hm AB for the WFC3/IR fil¬ 
ters. The {J\iQ-Hm) color can be shifted to the Vega system 
by adding 0.4919 mag. 

The results of our calibration analysis show that IR SBF 
measurements, especially in FllOW, can produce high- 
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accuracy distance measurements for red early-type galaxies 
with (g 475 -Z 85 o) > 1 -2 and (Jno-Hm) > 0.2 mag using Equa¬ 
tions ([^ through (j^. It is worth emphasizing again that the 
linear calibrations should not be extended to bluer colors: 
the bluer galaxies, comprising dwarf ellipticals and low-mass 
SOs, clearly do not follow extrapolations of the linear relations 
in Equations ([^ through (j^. Such blue galaxies likely have 
a wider variety of stellar populations than the giant ellipticals 
and SO galaxies, and thus exhibit larger scatter, even with re¬ 
spect to the quadratic fits; they are not well-suited for highly 
accurate distance measurements. The stellar population im¬ 
plications are discussed in detail below in Section]^ 

4.2. Independent Checks of the IR SBF Calibration 
Zero Point 

Now, as for several decades, the forefront of progress in the 
measurement of extragalactic distances is limited primarily by 
the uncertainty in the calibration zero point (e.g., Ereedman & 
Madore 2010; Riess et al. 2011). We have chosen to calibrate 
the WEC3 IR SBE distance scale using the extensive Virgo 
and Eornax optical SBE measurements made by Blakeslee et 
al. (2009) and their collaborators. This guarantees that the 
IR observational uncertainties and population variations will 
dominate the calibration uncertainty, not the precision of the 
reference distances. It does not, however, reduce the system¬ 
atic zero point uncertainty present in the optical SBF mea¬ 
surements, which in turn were based on HST Cepheid dis¬ 
tances (Freedman et al. 2001). Blakeslee et al. (2010) dis¬ 
cuss in detail the small offsets between several of the largest 
SBF surveys, including Tonry et al. (2001) and Jensen et al. 
(2003), and the application of metallicity corrections to the 
HST Cepheid distance scale of Freedman et al. (2001). The 
systematic uncertainty in the SBF distance scale due to the 
uncertainty in the Cepheid zero point is about 0.1 mag (Freed¬ 
man & Madore 2010; Blakeslee et al. 2010). 

One approach for avoiding the Cepheid zero point uncer¬ 
tainty would be to use theoretical stellar population model 
predictions to calibrate the absolute M in galaxies of varying 
ages, metallicities, and colors. This model-based approach 
would therefore make SBF a primary distance indicator inde¬ 
pendent of all other distance measurements, dependent only 
on our understanding of the luminosities and colors of red 
giants and other evolved stars of a particular age and metal¬ 
licity. Comparisons with several different stellar population 
models are presented below in Section]^ As will be shown, 
infrared stellar population models are not sufficiently consis¬ 
tent to provide a robust zero point for distance calibration at 
the 10% level. At present, we find that observed IR SBF 
magnitudes are more useful for constraining stellar population 
models than the models are for constraining the SBF distance 
calibration in the IR. 

Another approach is to find other distance indicators that 
are independent of the Cepheid calibration, such as the ge¬ 
ometrical distance to NGC 4258. Water masers orbiting the 
central black hole in NGC 4258 have now been used to ac¬ 
curately determine the distance to this galaxy using a purely 
geometrical technique based on the Keplerian orbits of the 
masers (Humphreys et al. 2013). While SBF magnitudes 
are best measured in early-type galaxies, archival //leo im¬ 
ages of the central bulge of NGC 4258 (GO-11570) provided 
us with an opportunity to explore the SBF calibration inde¬ 
pendent of the Cepheid distance scale. Given that an SBF 
measurement to NGC 4258 could allow us to bypass the sys¬ 
tematic uncertainty in the Cepheid calibration, we felt it was 


worth an attempt. Unfortunately, the presence of clumpy 
dust and recent star formation prevented us from achieving 
an accurate calibration using this galaxy, even when optical 
color images were used to identify dusty regions. The maser- 
calibrated fluctuation magnitude is significantly brighter (by 
~0.3 mag) than the calibration determined above for the el¬ 
liptical and SO galaxies in the Fornax and Virgo clusters with 
similar (g 475 -Z 85 o) colors (see Fig. 1^. This result was not a 
surprise; patchy dust adds to the fluctuation signal, as does 
the presence of younger populations containing bright AGB 
stars. Unfortunately, the geometrical distance to NGC 4258 
does not provide a useful direct calibration of the SBF tech¬ 
nique for elliptical galaxies. Given that optical Cepheid dis¬ 
tances to NGC 4258 have now been published by Macri et 
al. (2006), Fausnaugh et al. (2015), and Hoffmann & Macri 
(2015), the future value of NGC 4258 in calibrating IR SBF 
is therefore most likely to be through an improved calibration 
of the Cepheid distance scale zero point and metallicity cor¬ 
rections. 

Type la supernovae are one of the most accurate and 
widely-used distance measurement techniques in use today. 
In a recent paper, Cantiello et al. (2013) reported WFC3 /no 
and //i6o measurements of the SBF distance to NGC 1316, a 
type la supernova host galaxy in the Fornax cluster. They used 
the Jensen et al. (2003) F160W SBF calibration forNICMOS 
(including a metallicity correction to the Freedman et al. 2001 
Cepheid zero point) and applied a 0.2 mag offset to the NIC- 
MOS zero point to account for the difference in Alter width 
between NICMOS and WFC3/IR, based on predictions of the 
SPoT stellar population models (Raimondo 2009; Raimondo 
et al. 2005). To avoid uncertainties arising from differences 
between the methods used by Cantiello et al. and those used 
herein, we repeated the SBF analysis for NGC 1316 using the 
original WFC3/IR data (GO-11691) and the procedures de¬ 
scribed above. Our measured SBF magnitudes are listed in 
Tables and they are consistent with the Cantiello et al. 
(2013) values within the stated uncertainties. This galaxy is 
not an ideal SBF candidate due to the presence of extensive 
patchy dust near the center, but the regions farther out from 
the center appear relatively clean and the SBF signal is very 
strong. If we adopt the Blakeslee et al. (2009) z-band SBF 
distance (Table |^, the M values we And for NGC 1316 are 
brighter than the calibration prediction by ~0.3 to 0.5 mag. 
The Tonry et al. (2001) ground-based /-band SBF measure¬ 
ment of 31.66 ± 0.17 supports the Blakeslee et al. (2009) dis¬ 
tance for NGC 1316. If NGC 1316 is located at the same dis¬ 
tance as the core ellipticals in Fornax, then it appears to have 
a significant population of younger AGB stars, most likely the 
result of star formation that took place during a major merg¬ 
ing episode a few Gyr ago, that biases the IR SBF magni¬ 
tude. It may also have additional undetected dust contribut¬ 
ing to the fluctuations. On the other hand, the IR and optical 
SBF distances are inconsistent with the published type la su¬ 
pernova distances, which place NGC 1316 ~0.25 mag closer 
than the Fornax cluster core. If we adopt the la supernova 
distance modulus for the three normal supernovae published 
by Stritzinger et al. (2010) of 31.248 ±0.034±0.04 instead 
of the z-band SBF distance, our IR SBF measurements would 
be much more consistent with the elliptical galaxy calibration. 
Additional work is needed to reconcile the SBF and supernova 
distances to this galaxy. 

The IR SBF measurements of the maser host galaxy 
NGC 4258 and the supernova host galaxy NGC 1316 are not 
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('^110 -f^ieo) 

Fig. 4.— Fits to the absolute fluctuation magnitudes M as a function of (/iio“^i6o) color derived using individual distances (left) and average cluster distances 
(right). Red symbols indicate /fi 60 measurements (top set of points), and lower set of blue points are the 7jio measurements. Symbol definitions are the same as 
in Figure]^ 


TABLE 6 

Linear Calibration Coefficients 


Calibration 

a 

b 

X^/dof 

rms 

Distances 


Mno vs. (g-z) 

-2.946 ±0.015 

2.16±0.15 

0.48 

0.075 

dust 

11 

Ml 10 vs. (g-z) 

-2.935 ±0.017 

1.90 ±0.16 

0.86 

0.086 

dust 

12 

Ml 10 vs. (g-z) 

-2.961 ±0.032 

1.86 ±0.32 

1.23 

0.101 

indiv 

11 

Ml 10 vs. (g-z) 

-2.951 ±0.034 

1.61 ±0.34 

1.40 

0.105 

indiv 

12 

Mho vs. (J-H) 

-2.964 ±0.032 

6.7 ±0.9 

0.94 

0.092 

dust 

12 

Mho vs. (J-I/) 

-2.977 ± 0.047 

5.9±1.4 

1.76 

0.130 

indiv 

12 

Mi 60 vs. (g-z) 

-3.699 ±0.028 

2.13 ±0.27 

1.35 

0.114 

dust 

11 

M 160 vs. (g-z) 

-3.689 ±0.029 

1.90 ±0.28 

1.59 

0.114 

dust 

12 

M 160 vs. (g-z) 

-3.712 ±0.042 

1.76 ±0.42 

2.49 

0.138 

indiv 

11 

M 160 vs. (g-z) 

-3.704 ±0.043 

1.55 ±0.44 

2.45 

0.134 

indiv 

12 

M 160 vs. (J-H) 

-3.718±0.035 

7.1± 1.1 

1.35 

0.121 

dust 

12 

M 160 vs. (J-H) 

-3.731 ±0.050 

6.4±1.6 

2.44 

0.162 

indiv 

12 


Note. — Linear fits for (^ 475 —Z 85 o)> 1-2 or (7iio—^fj6o)>fi-22. For Vega magnitudes, 
subtract 0.7595 mag from ^\TC3/IR /no and 1.2514 mag from //leo AB. 


^ IC 1919 was excluded from the (^ 475 -^ 850 ) calibration because (g 475 -Z 85 o)< 1-2- The 
(/iiQ—7fi6o) color for IC 1919 is greater than 0.22, so it was included in the (/iiO“^i6o) 
calibration. We have included (^ 475 —Z 85 o) calibrations with and without IC 1919 to show 
the relatively small influence this one galaxy has on the (^ 475 —Z 85 o) calibration. The recom¬ 
mended calibration in Equations (3) and (4) exclude IC 1919. 


TABLE? 

Quadratic Calibration Coefficients 


Calibration 

a 

b 

c 

xVdof 

rms (mag) 

Distances 

Mho vs. (g-z) 

-2.991 ±0.004 

1.73 ±0.025 

4.55±0.12 

1.11 

0.115 

indiv 

Mho vs. (J-H) 

-2.977 ±0.011 

6.3 ±0.4 

59±6 

1.06 

0.129 

indiv 

Mi 60 vs. (g-z) 

-3.723 ±0.003 

1.54 ±0.022 

3.38 ±0.09 

2.17 

0.128 

indiv 

M,6o vs. (J-H) 

-3.733 ±0.030 

6.9 ±0.9 

59 ±15 

1.74 

0.149 

indiv 


Note. — Quadratic fits for all 16 galaxies. Use only when necessary for (g 475 -Z 85 o)< 1-2 or 

(•/iiO“^^i6o) <0-22. 
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sufficiently accurate at present to reduce the systematic un¬ 
certainty in the IR SBF distance scale calibration to less than 
the 10% that we currently inherit from the Cepheid calibration 
of the optical SBF distances used for this study (Blakeslee et 
al. 2009, 2010). These galaxies have significant patchy dust 
and recent star formation, both of which enhance the IR SBF 
signal over what is typically observed in quiescent elliptical 
galaxies. 

4.3. Comparison ofWFC3 and NICMOS SBF Magnitudes 

Eight galaxies in our current sample, mostly in the Fornax 
cluster, were included in the NICMOS NIC2 calibration of 
Mi6o (Jensen et al. 2003). A comparison of the SBF magni¬ 
tudes is shown in Figure!^ To make the comparison using AB 
magnitudes, we added li313 mag to the published NIC2 SBF 
magnitudes to shift from the Vega magnitude system used by 
Jensen et al. (2003) to the AB mag systemn 

The difference between NIC2 and WFC3/IR F160W SBF 
magnitudes shows a modest color dependence (left panel in 
Fig. 1^. We fitted the slope (the color term in the conversion 
of NlC2 SBF magnitudes to WFC3/IR) including uncertain¬ 
ties in both {Jiio-Hm) color and mieo SBF magnitudes; the 
rms for the fit is 0.055 mag and the per degree of freedom is 
0.66. The right panel of FigureJ^compares the NIC2 F160W 
apparent SBF magnitudes to WFC3/IR, with and without the 
color correction. The dashed line in the right panel of Fig¬ 
ure |^is not a fit, but shows the 45-degree perfect correlation 
line. SBF migo measurements made using NIC2 prior to the 
installation of the NICMOS cryocoolers may be compared to 
WFC3/IR measurements using the relation 

inm = 1M160.NIC2-1-67(71 io-J/i6o) +0.46 (7) 

(all AB magnitudes). Because we are comparing apparent 
fluctuation magnitudes directly, we do not need to be con¬ 
cerned about differences in distance scale calibrations be¬ 
tween the two studies or cameras. 

5. STELLAR POPULATION MODELS 

We turn now to how IR SBL measurements, with their sen¬ 
sitivity to red giant branch and intermediate-age AGB stars, 
can expose interesting differences between galaxies with dif¬ 
ferent star formation histories and better constrain single¬ 
burst stellar population models. 

The brightening of the SBL magnitudes in elliptical and SO 
galaxies at intermediate colors seen in the centers of Ligures]^ 
andffl and the subsequent drop in SBL brightness in the bluest 
dwarf ellipticals in this sample, provide powerful new con¬ 
straints for stellar population models, which have only been 
compared to redder galaxies in previous near-IR SBL studies. 
The scatter among the bluest galaxies in our sample is much 
larger than the observational uncertainties, and many of these 
galaxies exhibit significant radial gradients in IR fluctuation 
magnitude, color, or both. The current sample includes bluer 
and fainter galaxies than are typically targeted for SBL dis¬ 
tance measurements, and the breaking of the age-metallicity 
degeneracy in near-IR fluctuations provides a unique oppor¬ 
tunity to explore the stellar populations in these galaxies. 

Single-burst stellar population (often abbreviated SSP) 
models with constant age and metallicity are frequently used 
to interpret broad-band colors and other properties of unre¬ 
solved stellar populations in distant galaxies, including SBL 

http://www.stsci.edu/hst/nicmos/documents/handbooks/DataHandbookv8/ 

nic_ch5.9.3.htnil 


magnitudes. These models are usually calculated by integrat¬ 
ing collections of properly-weighted isochrones, and can be 
used to compute predicted SBL magnitudes directly without 
any need to link the apparent SBL magnitudes to an exter¬ 
nal distance scale calibration. Examples of theoretical SBL 
comparisons include Worthey (1993, 1994), Liu et al. (2000, 
2002), Blakeslee et al. (2001), Cantiello et al. (2003), Rai- 
mondo et al. (2005), Marm-Lranch & Aparicio (2006), Bis- 
cardi et al. (2008), and Lee et al. (2010). 

We compared our SBL measurements to three recent sets 
of SSP models for which 7iio and SBL magnitudes have 
been computed. The purpose of these comparisons is to ex¬ 
plore the limitations of our SBL distance calibration and pro¬ 
vide input to researchers working to improve stellar popu¬ 
lation models, particularly for understanding galaxy evolu¬ 
tion, when the observations may not so easily distinguish the 
effects of age and metallicity as our near-IR SBL measure¬ 
ments do. While we compared our IR SBL measurements to 
single-burst population models, real galaxies are composed 
of composite stellar populations with potentially many bursts 
of star formation. Because the fluctuations are dominated by 
the most luminous stars weighted as L^, they are even more 
strongly weighted towards young, luminous populations than 
are broad-band galaxy colors (Tonry & Schneider 1988). A 
composite population model will therefore predict an SBL 
magnitude close to that of the younger (or brighter) model 
component, even when the young population is only a small 
fraction (^10% to 20%) of the galaxy by mass (Jensen et al. 
2003; Blakeslee et al. 2001; Liu et al. 2002). The comparison 
to SSP model ages shown in this section should be considered 
the time since the most recent episode of star formation, not 
the average age of the dominant stellar population by mass. 

5.1. Teramo BaSTIModels 

The first set of stellar population models we consider here 
are based on the Teramo BaSTI model^ using a standard 
Salpeter initial mass function (IML) with ^ow-mass cutoff of 
0.5 Mq (Lee, Worthey, & Blakeslee 2010; Pietrinferni et al. 
2004, 2006; Cordier et al. 2007). We compared our observed 
fluctuation magnitudes Mno to absolute fluctuation magni¬ 
tudes computed for two variants of the BaSTI models: the 
solar-scaled abundance ratio models without convective over¬ 
shoot, and the a-enhanced version of the models described 
by Lee et al. (2010). The latter models have a mean [a/Le] 
of ~ 0.4 dex, but with physically-motivated variations in the 
abundance ratios of the individual a elements (i.e., O, Mg, Si, 
S, Ca, Ti); the [Le/H] abundances in these models have been 
correspondingly reduced to keep a fixed [Z/H]. 

The solar-scaled models shown in the top panel of Ligure|^ 
do not match the reddest giant ellipticals as well as the a- 
enhanced models shown in the lower panel (the solar-scaled 
models are included here to provide a point of reference 
with past SBL-model comparisons that only used solar-scaled 
metallicity models). Both model variations suggest that the 
intermediate-color galaxies in our sample have younger pop¬ 
ulations, as expected due to the presence of intermediate-age 
AGB stars. The models imply that the bluest galaxies are old 
and metal-poor. While we have chosen to show the compar¬ 
ison for Ml 10 and (g 475 -Z 85 o), the conclusions are the same 
when we compare the models to Mieo instead of Mno or 
(7 iio-Mi6o) instead of (g 475 - 285 o)- 

^ http://193.204.1.62/index.html 
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Fig . 5.— Comparison of NICMOS NIC2 and WFC3/IR /fi 60 SBF measurements. The left panel shows the fit to the color term used to correct the points in the 
right panel, which shows the corrected NIC2 SBF magnitudes plotted as a function of the WFC3/IR SBF magnitudes. For reference, the SBF magnitudes with 
no color con'ection are shown with X symbols. 



Fig. 6. — Teramo BaSTI models compared to Muo as a function of galaxy 
(? 475 “Z 85 o) color for two metallicity variants of the models computed using 
the individual galaxy distances: the “SSS” models (top panel) are solar-scaled 
abundances spanning the range from [Fe/Ft]= —0.659 to -1-0.395. The “AES” 
models (lower panel) have enhanced a element compositions and reduced 
[Fe/Fl] to make the overall metallicity Z the same as the solar-scaled models 
shown in the top panel. The range of [Fe/H] spanned by the AES models 
is —1.01 to -1-0.05. Lines of constant metallicity are shown with red dashed 
lines, and lines of constant age from 2 to 13 Gyr are plotted with dotted blue 
lines. 

5.2. Teramo SPoTModels 

The second set of models was developed by the Teramo 
SPoT group (version BaSeLS.l, Raimondo 2009; Raimondo 
et al. 2005). The Teramo SPoT models pay special attention 
to the thermally-pulsating asymptotic giant branch (TP-AGB) 
stars, particularly in young to intermediate-age populations. 
The SPoT model SBF magnitude predictions in the near-lR 
have been empirically compared to a variety of clusters in the 
Large Magellanic Cloud and have been shown to match em¬ 
pirical measurements of SBF magnitudes, integrated magni¬ 
tudes, star counts, and colors (Raimondo 2009; Cantiello et 
al. 2007; Cantiello 2012). 

The SPoT models in Figure 1^ differ significantly from 
the BaSTI models in Figure ^though they are based on 
the same stellar evolutionary tracks as the BaSTI models in 
the upper panel of Fig. The SPoT models derive SBF 
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Eig. 7.— Recent versions of the Teramo SPoT models compared to Mho 
measurements as in Figure!^ The metallicity range for these models spans 
[Fe/Fl] = —0.66 to -1-0.4 and the ages from 2 to 14 Gyr. The lines and symbols 
are defined in the same way as in Figures|^and|^ 

magnitudes in a procedure that allows modelers to statisti¬ 
cally combine various stellar population models produced by 
stochastic variations in the number and properties of bright 
and rare stars, including the TP-AGB and horizontal branch 
populations in intermediate-age and old stellar populations. 
Overall, the SPoT and BaSTI models agree for the red and 
intermediate-color ellipticals, with somewhat fainter SBF 
magnitudes for the oldest and most metal-rich giant ellipti¬ 
cals. The SPoT models predict a larger spread in SBF mag¬ 
nitude at younger ages and imply that the bluer galaxies in 
the sample are all younger than about 5 Gyr, in contrast with 
the BaSTI models, which span a narrower range in M\\q and 
imply ages greater than 5 Gyr for the bluest galaxies. 

5.3. Padova Models 

The third set of models are based on the Padova isochrones 
(Fig .[^, which include sophisticated handling of the TP-AGB 
evolutionary phase. The Padova tracks use solar-scaled metal¬ 
licity abundance ratios and do not include a-element enhance¬ 
ment. SBF magnitudes were computed by Lee et al. (2010) 
using the evolutionary tracks of Marigo et al. (2008). 

The Padova SSP models show the largest spread in Mno 
of all the models considered here. At the red end, the giant 
ellipticals again agree with old, metal-rich population mod- 
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Fig. 8 .— Stellar population models based on the Padova isochrones com¬ 
pared to Ml 10 measurements computed usingindividual SBF distances. The 
symbol definitions are the same as in Figured Lines of constant metallicity 
from [Fe/H] = - 0.705 (in the extreme lowerleft comer) to 0.0 and +0.222 
(overlapping at right) ai'e shown with red dashed lines. Lines of constant age 
from 2 to 13 Gyr ai'e plotted with dotted blue lines. 

els, with the intermediate bluer galaxies having brighter fluc¬ 
tuations and younger populations. The bluest galaxies in 
the sample are consistent with ages intermediate between the 
BaSTI and SPoT models, in the range 3 to 7 Gyr. Because of 
their larger spread in SBF magnitude predictions, the Padova 
models, as compared to the other sets of models, suggest 
that the data provide better discrimination between ages and 
metallicities. 

5.4. Fluctuation Colors 

SBF measurements at two wavelengths can be used to elim¬ 
inate uncertainties resulting from distance error; for instance, 
the “fluctuation color” (1 miio-1mi6o) can be compared to stel¬ 
lar population models in a distance-independent way. Figure!^ 
compares ((Mno-iwieo) to the Teramo BaSTI solar-scaled and 
a-enhanced models as a function of (g 475 -Z 85 o)- As discussed 
by Lee et al. (2010), the predicted SBF magnitudes are sen¬ 
sitive to a-element abundance mostly because of the effects 
of oxygen-enhancement on the upper red giant branch and 
AGB phase. The data in Figure 1^ agree on average with the 
a-enhanced models relatively well, but only poorly with the 
solar-scaled models. However, this set of a-enhanced mod¬ 
els also predicts a narrower range of (miio-mieo) than is ob¬ 
served. In particular, the low-mass dwarf IC 3032 in Virgo 
agrees better with the locus of the solar-scaled models; this 
may indicate real variation in a-element abundance ratios 
among the sample galaxies. We conclude that fluctuation col¬ 
ors can provide useful information on elemental abundance 
ratio trends with age and metallicity in elliptical galaxies, in¬ 
dependent of the uncertainties in the distance calibration. In 
addition, comparisons with other observables that are sen¬ 
sitive to age, metallicity, or a-enhancement should provide 
powerful joint constraints for future stellar population mod¬ 
els. 

5.5. Radial SBF Gradients 

To further explore the origins of the scatter in SBF mag¬ 
nitude among the low-luminosity blue dwarf ellipticals, we 
measured the radial behavior of the SBF amplitude. Many 
of the galaxies in our sample, particularly the low-luminosity 
galaxies, are quite elliptical. To get the cleanest gradient mea¬ 
surements possible, we repeated the SBF analysis and mea¬ 
sured (g 475 -Z 85 o) colors in elliptical annuli for a subset of the 
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Fig. 9.— Distance-independent IR fluctuation colors for the solar-scaled 
and a-enhanced Teramo models. Blue dotted lines indicate lines of con¬ 
stant age, from 2 to 13 Gyr. Red dashed lines show constant metallicity 
tracks, from [Fe/Fl] = —0.659 at the left to +0.395 at right (upper panel) and 
[Fe/H] = — 1.01 to +0.05 (lower panel). The two points that are not labeled for 
clarity are NGC 1399 in Fornax and NGC 4636 in Virgo. Symbol definitions 
are the same as in Figures|^and[^ 



Fig. 10.— Radial SBF and color gradients for the survey galaxies, with 
larger symbols corresponding to the inner annuli in which SBFs were mea¬ 
sured and smaller symbols being outer regions. Black lines link measure¬ 
ments in the same galaxy. The measurements are compared to the BaSTI 
a-enhanced models; the conclusions are the same for the other models be¬ 
cause they all have lines of constant metallicity with similar slopes for this 
color range. 1C 2006 is not labelled but falls very close to NGC 1374. 

galaxies with significant M gradients (IC 1919, NGC 1375, 
NGC 1380, IC 3025, IC 3487, and IC 3586; see Fig.|^ For 
the rest of the sample we used the circular annulus SBF and 
(g 475 -Z 85 o) color measurements previously used for the cali¬ 
bration. The results are plotted in Figure [10| 

The comparison of radial SBF gradients to the models 
shows a significant difference between the giant ellipticals 
and the smaller galaxies. Most of the galaxies show a gra¬ 
dient sloping from upper left to lower right in Figure [T^ 
roughly along lines of constant metallicity. While the vari¬ 
ous models look quite different in detail, they all have lines 
of constant metallicity sloping in approximately the same di¬ 
rection at the relevant colors. The lower-luminosity galax¬ 
ies have centers (larger symbols in Fig. consistent with 
younger stellar populations (brighter fluctuations) and nearly 
constant metallicities: IC 1919, IC 3025, IC 3487, IC 3586, 
and NGC 1375 all have brighter fluctuations near their cen- 
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ters. In contrast, the giant ellipticals on the red end of the 
diagram tend to show older, and sometimes more metal-rich, 
populations in their centers: IC 3032, IC 2006, NGC 1374, 
NGC 1399, NGC 1404, NGC 4458, and NGC 4649 have 
fainter fluctuations towards their centers. Four of the sixteen 
galaxies (NGC 1380, NGC 1344, NGC 4458, and NGC 4472) 
appear to have color and SBF gradients consistent with little 
or no age variation. The majority of the galaxies show gradi¬ 
ents consistent with primarily age variations, however. 

The IR SBF and color gradients hint at differing formation 
histories for the galaxies. As noted above, low-luminosity 
blue elliptical galaxies usually have older populations at large 
radii, and thus to have formed stars more recently near their 
centers from metal-poor gas, while the giant ellipticals formed 
stars in their cores long ago from enriched gas. Optical studies 
of SBF gradients, in contrast, suggest that the outer parts of 
giant ellipticals have colors and fluctuation magnitudes con¬ 
sistent with lower metallicity populations and little or no age 
gradient (Cantiello et al. 2005, 2007). Our IR measurements 
do not cover as large a range in radius as the optical studies 
because of the smaller held of view spanned by the IR detec¬ 
tors and the brighter sky background. Because the data here 
are confined within the effective radius of the several largest 
galaxies (Table [T]), we cannot derive strong constraints on the 
large-scale stellS population gradients in these giant galax¬ 
ies. Interestingly, the physical size (r <5 kpc) of the radial 
region probed by our SBF data is similar to that of the lumi¬ 
nous dense cores observed at high redshift (e.g., van Dokkum 
et al. 2010; van der Wei et al. 2014), which are believed to be 
the seeds around which massive modern-day ellipticals were 
assembled through hierarchical merging and accretion. Our 
results on the old ages and relative homogeneity of the stellar 
populations of the giant ellipticals within the central region 
are thus in line with this scenario for early-type galaxy evo¬ 
lution, if the elliptical galaxies in our sample may be com¬ 
pared with the high-redshift counterparts (Fritz et al. 2014). 
At much larger radii in giant ellipticals, evidence exists from 
optical and IR color gradients that the stellar populations at 
several effective radii are older than the mean stellar age in 
the core (La Barbera et al. 2012; but see also Greene et al. 
2015). While our data do not probe such large radii in the 
massive galaxies, we And the same trend for the lower mass 
galaxies where our measurements do reach beyond the effec¬ 
tive radius. It would be interesting to extend the IR SBF gra¬ 
dients to correspondingly large radii in the massive ellipticals 
to see if these also show evidence for larger ages in the galaxy 
outskirts. 

Crucially for the purpose of distance measurement, the ra¬ 
dial and population variations in the reddest galaxies cause 
SBF magnitude variations along the direction of the linear cal¬ 
ibration. These galaxies are ideally suited for distance mea¬ 
surement because variations from galaxy to galaxy in age or 
metallicity are adequately accommodated by the calibration 
slope with (g 475 -Z 85 o)- On the other hand, the radial gradi¬ 
ents in SBF magnitude in the bluer dwarf galaxies (the gray 
symbols in the left half of Fig. [TOl are perpendicular to the 
quadratic calibration relation, leading to greater scatter in the 
distance calibration. These low-luminosity dwarfs show a 
wide range of ages and metallicities, with several showing ev¬ 
idence of recent star formation near their centers, as revealed 
by their SBF magnitudes. 

Only one set of SBF measurements and models is shown 
in Figure [T0{ the general conclusions, however, are consistent 
for all the models discussed herein. The radial variations are 



Fig. 11.— Lick index measurements for eleven of the galaxies in our sam¬ 
ple (Kuntschner 2000; Trager et al. 2000; Caldwell et al. 2003) as compared 
with the ck-enhanced models from Lee & Worthey (2005). H/3 is more age- 
sensitive, while y[gb is metallicity-sensitive. Such data provide an indepen¬ 
dent way of checking the SBF model comparison conclusions; see text. 


consistently sloped along lines of constant metallicity for all 
the models, so even though the models might not agree on the 
specific age and metallicity of a particular galaxy, the trend 
towards younger ages in the centers of bluer dwarf galaxies is 
consistent for all the models. 

5.6. Line Index Age and Metallicity Constraints 

Eleven of the galaxies in our sample have published HP 
and Mgb Lick index measurements from Kuntschner (2000), 
Trager et al. (2000), and Caldwell et al. (2003). The former 
index is more sensitive to age, and the latter to metallicity. 
Figure [^compares these data with the -1-0.3 dex a-enhanced 
model predictions from Lee & Worthey (2005). The line in¬ 
dex measurements generally sample different regions of the 
galaxies than our SBF data, but we see again that the red¬ 
dest ellipticals agree with the old, metal-rich population mod¬ 
els. The youngest galaxies as determined using absorption 
lines do not always agree with the SBF models. IC 3487 
has the strongest H/3 index but the faintest Muo among the 
bluer galaxies, implying a relatively older age compared to 
NGC 1375, which has somewhat weaker H/3 and brighter 
Ml 10 (although there is considerable variation between SBF 
models at the youngest ages). Further work is needed to ex¬ 
plore the radial behavior of Mgb and H/3 as compared to our 
SBF measurements on the same scales, particularly since we 
detect a significant radial age gradient in many of the lower- 
luminosity galaxies. 

5.7. Model Comparison Summary 

The SBF models shown in Figures [^|7] and|^agree for old, 
metal-rich populations such as those commonly found in gi¬ 
ant elliptical galaxies, and for which extensive comparisons 
have been made in the past (e.g., Jensen et al. 2003; Cantiello 
et al. 2012; Fritz 2012). At younger ages and lower metallici¬ 
ties, the three sets of models are significantly different, and 
the conclusions we draw about the ages of the blue dwarf 
ellipticals are strikingly different. These three sets of mod¬ 
els allow for the blue, low-mass dwarf ellipticals to poten¬ 
tially have a range of ages from 2 to 14 Gyr, and [Fe/H] from 
about -1 to -0.3. Optical /-band SBF measurements of blue 
dwarf ellipticals (Mieske et al. 2006) also show a large spread 
of fluctuation magnitudes and younger implied ages. Since 
each set of models was computed with the aim of understand- 
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ing something different—the role of a-enhancement in early- 
type galaxies for the BaSTI models (Lee et al. 2010) and the 
TP-AGB phase for the SPoT and Padova models (Raimondo 
2009; Cantiello et al. 2007, 2012; Lee et al. 2010)—the differ¬ 
ences shown here provide the starting point for future detailed 
comparison of the near-IR properties of actual red giant and 
AGB stars in unresolved stellar populations with model pre¬ 
dictions. We have chosen to show model comparisons using 
Ml 10 plotted against (g 475 -Z 85 o) computed using the individ¬ 
ual galaxy optical SBF distances. The general trends and con¬ 
clusions are similar when the models are compared to Mieo or 
plotted against (Jno-Him). 

It is beyond the scope of this study to provide a critical anal¬ 
ysis of the strengths and shortcomings of each of these sets of 
models. For the purposes of this study, we conclude that the 
IR SBF distance calibration is robust when applied to old, red, 
metal-rich galaxies like the giant ellipticals typically found in 
environments like Virgo and Fornax, and even to some with 
intermediate-age populations and somewhat bluer colors. The 
bluer dwarf ellipticals, on the other hand, provide important 
new observational constraints that should be of interest to re¬ 
searchers constructing the next generation of stellar popula¬ 
tion models. It is clear that better constraints from data such 
as these will be valuable in better defining the properties of 
young, metal-poor populations, and the brightness of the TP- 
AGB stars within these populations. 

6. RECOMMENDATIONS EOR MEASURING IR SBE 
DISTANCES WITH WEC3 

WEC3 on the HST makes it possible to measure IR SBF 
magnitudes at relatively large distances in modest exposure 
times. Based on our experience with the Fornax and Virgo 
cluster calibration data presented here, as well tests with 
WFC3/IR observations in the Coma cluster from HST pro¬ 
gram GO-11711 (see Blakeslee 2013) and with the instrument 
exposure time calculator, we provide the following guidelines 
to help other astronomers plan WFC3 SBF observations and 
make use of existing data in the HST archive. 

1. The fluctuations are brighter at //leo than at 7no, but the 
7iio filter is signficantly wider than Hi(,o (by 0.8 mag), 
which largely cancels out the brightness advantage. The 
image quality is slightly better at Ji lo than //igo, and the 
background is slightly lower. The net effect is that ex¬ 
posure times to achieve a particular SBF S/N ratio is 
about the same in the two filters, but the ability to de¬ 
tect and remove contaminating point sources (primar¬ 
ily globular clusters) is better at /no, and there is less 
scatter in the calibration. We therefore advise choos¬ 
ing JiiQ over //i6o when possible. The broad F140W 
filter would be a good alternative, unhampered by the 
1.083 fim He emission line in the upper atmosphere, but 
that bandpass remains uncalibrated for SBF. 

2. Typical one-orbit exposure times (~2400 s) suffice for 
measuring distances in either /no or //igo out to about 
80 Mpc. This distance limit is imposed by the point 
source sensitivity required to detect and remove glob¬ 
ular clusters from the image, with the goal of reaching 
within ~ 1 mag of the peak of the globular cluster lu¬ 
minosity function in /no, or ~ 1.5 mag of the peak in 
//i6o, for which the fluctuations are relatively brighter 
(Jensen et al. 1998). The SBF signal itself can be de¬ 
tected to much larger distances (beyond 100 Mpc) in a 


single orbit, but the large correction for contaminating 
point sources would then dominate the uncertainty. 

3. Exposure times for more distant galaxies should be 
scaled to achieve a point source sensitivity sufficient 
to detect and remove globular clusters 1 to 1.5 mag 
brighter than the peak of the globular cluster luminos¬ 
ity function. The exposure time needed to detect the 
stochastic fluctuations that comprise the SBE signal 
goes as d^, but the time required to detect the globu¬ 
lar clusters increases significantly faster, scaling as 

to d‘^, because of the bright background on which the 
globular clusters are superimposed. 

4. To avoid issues with correlated noise, do not use the 
defaultpipeline-combined images. Use the^f files 
without correcting the WEC3/IR field distortion. It may 
be possible to recreate the drz files using the lanczos3 
kernel in the astrodrizzle package; this approach was 
not tested for WEC3/IR as part of this study, but has 
been used successfully in the past for ACS data (cf. 
Cantiello et al. 2005; Mei et al. 2005a). 

5. Because SBE magnitudes depend on the properties of 
the stellar populations, high-quality color data are es¬ 
sential to determine accurate distances. Eor the most ac¬ 
curate SBE distances, one should target giant elliptical 
and SO galaxies with old stellar populations, for which 
(^ 475 - 2 s 5 o) colors are greater than 1.2 or (/no-Migo) are 
greater than 0.22 AB mag. The population variations 
at bluer colors are too great for robust distance mea¬ 
surements. If possible, it is best to use ACS (g 475 -Z 85 o) 
colors, but WEC3/IR {Ji io-Hi(,o) colors are an accept¬ 
able alternative. If necessary, other color indices may 
be translated to (^ 475 -^ 850 ) or {Jiw-Hm) using well- 
constrained empirical or model relations for old, metal- 
rich populations. 

7. SUMMARY 

We have measured Jno and //leo SBE magnitudes and 
(7 iio-//i6o) colors for 16 early-type galaxies in the Virgo and 
Eornax clusters observed with WEC3/IR. All of these galaxies 
have previously measured SBE distances and (g 475 -Z 85 o) col¬ 
ors from ACS. We find that the luminous red galaxies in the 
sample follow linear relations between absolute SBE magni¬ 
tude and optical or IR color; SBE distances to such galaxies 
can be determined within a statistical uncertainty of 5% us¬ 
ing the calibration relations that we have presented in Equa¬ 
tions ([^ through (j^. The systematic uncertainty of this cali¬ 
bration IS ^ 10% due to the uncertainty in the Cepheid calibra¬ 
tion on which this work is based. Stellar population models 
are not consistent enough to provide a direct calibration of the 
IR SBE technique accurate to 10%, but they are valuable for 
inferring age and metallicity trends in and among the sample 
galaxies. The scatter in SBE magnitude among bluer galax¬ 
ies is large, indicating that a wider variety of stellar popula¬ 
tions dominate the light from these galaxies and that a simple 
broadband color does not adequately parameterize the com¬ 
plexities inherent in such populations. Some galaxies, partic¬ 
ularly those of intermediate color, clearly have younger (4 to 
8 Gyr) populations, likely with AGB stars that enhance the IR 
fluctuation amplitude. Bluer galaxies, primarily dwarf ellip¬ 
ticals, may have very-low metallicities and/or younger ages, 
with the interpretations varying among different sets of stellar 
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population models. Finally, we have provided practical advice 
to guide researchers interested in undertaking SBF measure¬ 
ments with WFC3/IR. 
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